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ABSTRACT 


This  study  demonstrates  the  significance  of  water  vapor's 
influence  on  satellite-retrieved  aerosol  characteristics  using 
NOAA's  AVHRR.  An  improvement  to  optical  depth  and  Aerosol 
Particle  Size  Index  (Si 2)  estimations  derived  from  channels  1 
(0.63pm)  and  2  (0.86pm),  is  made  through  knowledge  of  column 
water  vapor  derived  from  channels  4  (10.8pm)  and  5  (12.0pm). 
A  2.0  gm  cm'2  column  water  vapor  produces  a  5%  increase  in 
S12.  This  results  in  a  15%  error  in  the  variable  scattering  phase 
function  P(0)  and  retrieved  aerosol  optical  depth  (5a).  The  error 
introduced  by  water  vapor  is  quantified  through  use  of  the 
LOWTRAN7  atmospheric  propagation  model  to  be  applied  as  a 
reformulated  parameterization  of  (P(0)). 
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I.  INTRODUCTION 


A.  PURPOSE 

Direct  measurements  of  aerosol  characteristics  can  be  provided  by  ship  or 
aircraft.  Although  they  provide  valuable  information,  they  are  very  expensive, 
require  detailed  coordination,  and  only  small  sections  of  the  world's  oceans  can 
be  monitored  at  any  time.  Satellites  can  provide  near  continuous  spatial  and 
temporal  data  desirable  for  evaluating  large  scale  global  effects  as  well  as  timely 
local  measurements  suitable  for  military  applications.  Meteorological  satellites, 
however,  have  been  designed  to  measure  cloud  albedos  while  much  lower 
radiances  are  associated  with  aerosol  optical  depth  variations  (Durkee  1984). 
Therefore,  careful  processing  of  satellite  data  is  necessary  to  determine  aerosol 
characteristics. 

Frost  (1988)  has  developed  an  algorithm  to  derive  aerosol  particle 
characteristics  from  the  Advanced  Very  High  Resolution  Radiometer  (AVHRR) 
on  the  NO  A  A  polar  orbiting  satellites.  By  taking  advantage  of  the  varying 
sensitivities  of  aerosol  scattering  in  two  spectral  windows,  a  multi-channel 
comparison  provides  information  about  the  aerosol  particle  size  distribution.  By 
applying  this  information  to  a  variable  scattering  phase  function,  improved 
estimations  of  optical  depth  are  possible. 

This  study  will  provide  an  improvement  to  Frost’s  technique.  Frost  requires 
analyzed  pixels  to  be  free  of  clouds.  However,  these  clear  pixels  still  contain 
varying  amounts  of  water  vapor.  Frost  assumes  that  varying  amounts  of  water 
vapor  do  not  appreciably  effect  retrieved  radiances.  This  study  will  use  the  U.  S. 
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Air  Force's  LOWTRAN7  atmospheric  propagation  model  to  demonstrate  the 
degree  that  water  vapor  influences  present  aerosol  estimations. 

This  thesis  will: 

1 .  Show  the  effect  of  varying  water  vapor  amounts  on  satellite- 
retrieved  aerosol  characteristics,  specifically  the  particle  size  index 
parameter  or  S12  as  defined  by  Frost  (1988). 

2.  Provide  a  parameterization  of  S12  corrected  for  column 
water  vapor  (w),  using  a  multichannel  analysis  of  AVHRR  data. 

B.  MOTIVATION 

Radiation  in  the  atmosphere  is  attenuated  by  many  processes.  It  is  absorbed 
and  scattered  by  molecules,  aerosols,  and  meteorological  features  such  as  clouds, 
fog,  and  precipitation.  Some  of  these  influences  may  be  of  short  duration  and  of 
local  scale  such  as  a  cumulus  cloud  passing  overhead  and  casting  a  shadow. 
Others,  such  as  stratospheric  dust  particles  are  semipermanent  and  of  global 
proportions. 

Aerosols  are  solid  or  liquid  particles  dispersed  in  the  atmosphere.  Dust, 
smoke,  and  even  fog  or  rain  are  examples  of  aerosols.  Aerosol  particles  are 
present  in  the  atmosphere  in  various  forms  and  concentrations.  In  this  thesis,  the 
change  in  particle  size  due  to  coalescence  of  water  vapor  will  not  be  considered. 

With  the  ever  increasing  use  of  electro-optical  devices,  military  applications 
benefit  significantly  from  improved  aerosol  measurements.  The  optimum 
waveband  for  a  laser  is  from  0.4  to  1.0  pm.  This  waveband  avoids  ozone 
absorption  at  shorter  wavelengths,  as  well  as  the  water  vapor  absorption  that 
dominates  longer  wavelengths  (Bloembergen,  et  al  1987).  This  waveband, 
however,  is  significantly  influenced  by  aerosol  scattering. 
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Another  important  hypothesis  that  demonstrates  the  need  for  better  aerosol 
estimations  and  modeling  on  a  much  larger  scale  is  that  of  Charlson,  et  al. 
(1987).  Charlson  proposes  that  the  major  source  of  cloud  condensation  nuclei 
(CCN)  over  remote  unpolluted  oceans  appears  to  be  aqueous  dimethylsulphide 
gas  (DMS)  produced  by  phytoplankton  as  a  waste  product.  DMS  is  emitted 
through  the  sea  surface  through  a  chain  of  complex  biological  and  chemical 
processes  where  it  is  oxidized  in  the  atmosphere  to  form  sulfate  aerosol  particles 
which  act  as  CCN.  Increased  phytoplankton  populations  and  the  resulting 
increase  in  aerosol  particles  effectively  cool  the  atmosphere  by  increasing 
backscatter  to  space.  A  30%  increase  in  aerosols  over  the  ocean  surfaces 
worldwide  can  conservatively  decrease  the  average  global  temperature  by  1.3  K 
(Charlson  1987). 

For  many  applications,  aerosols  result  in  noise  interfering  with 
measurements  of  other  processes.  These  applications  benefit  from  improved 
aerosol  estimations  that  eliminate  aerosol  contamination  of  the  data. 

C.  THE  MODEL 

LOWTRAN7  (Low  Resolution  Transmittance)  is  a  model  that  calculates 
atmospheric  transmittance,  scattering  and  absorption.  It  incorporates  separate 
molecular  profiles  for  all  major  as  well  as  13  minor  and  trace  gases,  all  as  a 
function  of  altitude  (0  to  100km).  Representative  atmospheric  aerosol,  cloud, 
and  rain  models  are  provided  with  the  ability  to  replace  any  or  all  values  with 
theoretical  or  measured  data  at  the  discretion  of  the  user.  The  model  also 
provides  six  reference  atmospheres  with  options  to  redefine  or  retain  profiles. 

LOWTRAN7  is  a  band  model.  Average  transmittance  is  calculated  over  20 
cm-1  bands.  The  band  parameters  have  been  determined  semi-empirically 
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determined  semi-empirically  through  the  averaging  of  line  parameters.  In 
principle,  the  line  by  line  method  is  available,  but  the  calculations  are  too 
elaborate  and  time  consuming  for  most  users.  Some  accuracy  is  gained  by  a  line 
method,  but  it  is  not  significant  enough  to  warrant  the  loss  of  efficiency. 
Individual  lines  within  a  band  may  vary  significantly,  yet  it  is  fair  to  characterize 
the  effective  band  parameters  through  the  averaging  of  line  parameters  so  long 
as  a  specific  line  result  is  not  desired. 

For  use  in  multiwavelength  applications  (ie.  satellites)  the  appropriate 
spectral  response  functions  for  the  sensor  can  be  applied  to  the  output  data 
through  LOWFIL7  (LOWTRAN7's  filter  program).  The  data  set  derived 
through  the  running  of  the  main  program  is  filtered  through  the  sensor's 
response  function  and  produces  a  normalized  radiance  estimate  of  the  sensor 
measurement. 
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II.  THEORY 


A.  PARTICLE  SIZE  INDEX  (S12) 

The  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  has  been  a  part 
of  all  the  NOAA  polar  orbiting  satellites.  It  has  five  channels: 


TABLE  1 


NOAA-7  AVHRR  CHANNEL  BANDWIDTHS 


Channel  # 
1 
2 

3 

4 

5 


Bandwidth  (pml 
0.58  -  0.68 
0.725  -  1.10 
3.55  -  3.93 
10.30-11.30 
11.50-12.50 


AVHRR  satellite  detected  radiances  in  the  red-visible  (chan.  1)  and  the  near 
infra-red  (chan.  2)  wavelength  ranges  have  been  found  to  be  positively 
correlated  with  aerosol  optical  depth  (8)  (Durkee  1984;  Griggs  1979). 

Using  a  single  scattering  approximation  and  knowing  that  optical  depths  over 
the  ocean  are  small  (Pfeil  1986),  Frost  (1988)  used  a  simplified  form  of  the 
Radiative  Transfer  Equation  (RTE): 


L(8,M>)  -  — p(0)  8  (2.1) 

where,  L  =  diffuse  radiance, 

coo  =  single  scatter  albedo. 
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F0  =  solar  radiative  flux, 

(i  =  cos0  (0  =  satellite  zenith  angle), 
p  =  scattering  phase  function, 

0  =  single  scattering  angle, 

8  =  optical  depth. 

In  an  optically  thin  (small  8)  atmosphere,  multiple  scattering  effects  are  very 
small,  making  single  scattering  a  reasonable  assumption  (Frost  1988). 

There  are  three  primary  sources  of  upwelled  radiance  (L)  that  contribute  to 
satellite  detected  radiance  in  the  red-visible  and  near-infrared  wavelengths.  The 
contributions  can  be  expressed  as: 

L  =  La  +  Lr  +  Ls  (2.2) 

where  La  is  aerosol  scatter,  Lr  is  molecular  or  Rayleigh  scatter,  the  wavelength 
of  radiation  is  much  greater  than  the  size  of  the  scattering  material,  and  Ls  is  sea 
surface  reflectance.  Except  for  sunglint,  Ls  is  very  small  for  red  wavelengths 
(albedo  =  0.5%)  and  zero  for  wavelengths  greater  than  0.7|j.m  (Ramsey  1968). 
Therefore,  as  long  as  sunglint  geometry  is  avoided,  total  radiance  is  due  to  the 
sum  of  aerosol  and  molecular  scattering  effects  (L  =  La  +  Lr). 

Lr  is  strongly  wavelength  (A,)  dependant,  but  does  not  vary  spatially. 
Therefore  Lr  forms  a  baseline  L(X)  and  aerosol  particles  are  the  dominant 
source  of  upwelled  radiance  variations  (AL  ~  ALa).  From  Eqn.  2.1,  changes 
in  satellite  measured  radiance  are  directly  determined  from  changes  in  aerosol 
characteristics  like  coo,  p(0),  and  aerosol  optical  depth  (8a).  Therefore,  to 
retrieve  a  more  accurate  8  from  La  measurements,  knowledge  of  p(0) 
variations  is  required. 
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Durkee  (1984)  showed  that  retrieved  radiances  at  red  wavelengths  (chan.  1) 
are  more  sensitive  to  varying  densities  of  small  aerosol  particles  than  that  of 
near-infrared  (nir)  wavelengths  (chan.  2).  The  most  significant  parameter 
contributing  to  retrieved  radiances  at  AVHRR  channel  1  and  2  wavelengths  is 
the  extinction  coefficient  ((Text)  and  specifically  the  vertical  integral  of  cext  or 
the  optical  depth  (6). 


H 

8=d<^extdz  (2.3) 

0ext  is  the  sum  of  extinction  due  to  absorption  (Gabs)  and  due  to  scattering 
(tfscat)-  For  marine  aerosols,  absorption  is  negligible,  therefore  oext  =  Gscat 
and 


H 

5=  J<?scat  dz.  (2-4) 

Three  terms  contribute  to  aScat: 


Oscat  =  J  Jtr2  Qscat(m,X,r)  n(r)  dr  (2.5) 

where,  itr2  =  particle  cross  sectional  area, 

Qscat(mA^)  =  scattering  efficiency, 
m  =  complex  index  of  refraction, 

X  =  wavelength, 
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r  =  particle  radius, 


n(r)  =  number  distribution  of  particles. 

Figure  1  shows  the  influence  of  the  terms  of  Eqn.  2.5  as  a  function  of  r  and 
the  resulting  number-weighted  extinction  (7tr2Qscatn(r))  for  channel  1  and 
channel  2.  Scattering  cross  section  (rcr2  Qscat)  and  the  particle  size  distribution 
(n(r))  for  a  modeled  distributic of  marine  particles  at  80%  RH  are  plotted.  As 
wavelength  increases,  so  does  the  r  of  most  efficient  scattering  (Qscat)-  The 
result  is  a  shift  in  the  scattering  cross  section  and  a  similar  shift  in  the  r  that  most 
significantly  contributes  to  cScat.  and  therefore  L  at  a  given  k.  Therefore,  the 
slope  of  L(X.)  closely  resembles  the  slope  of  5(A.),  and  is  similar  to  the  slope  of 
the  aerosol  particle  size  distribution  (n(r)).  Since  aerosol  particles  most 
effectively  scatter  radiation  at  a  wavelength  similar  to  the  radius  of  the  particle, 
an  increase  in  the  number  of  small  aerosol  particles  will  increase  the  ratio  of 
LChl  to  Lch2.  and  the  spectral  variation  can  be  quantified: 


L(red) 

L(nir) 


L(chl) 

L(ch2) 


CQq’Fq 


p(0)-8(chl) 


^P(e>-s(ch2) 


5(red) 

5(nir) 


(2.6) 


This  ratio  is  called  the  Particle  Size  Index  (S12)  (Frost  1988):1 


„  _  (L)red  _  (S)red 

^  (L)nir  (5)nir 


(2.7) 


JWhile  Frost  removed  Lr  prior  to  taking  the  ratio,  it  will  not  be  removed  in 
this  study.  It  shall  be  accounted  for  in  the  parameterization  of  p(0). 
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Fig.  1  Extinction  formulation  as  a  function  of  particle  radius  for 
channels  1  (solid)  and  2  (dash).  Modified  from  Durkee  (1984). 
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Satellite  measured  radiances  therefore  cannot  provide  the  actual  particle  size 
distribution,  but  the  S12  gives  the  slope  of  the  size  distribution  curve.  In 
general,  the  greater  the  contribution  of  smaller  particles,  the  larger  the  ratio. 

Figure  2  shows  a  plot  of  NOAA-7  AVHRR  channel  1  and  channel  2 
response  function  curves  superimposed  on  LOWTRAN7's  Tropical  model 
atmosphere  output  radiances  (L(X)).  Although  as  this  paper  will  show,  water 
vapor  absorption  does  effect  S12,  the  dominant  variable  influencing  L(X)  in 
these  wavelength  bands  is  aerosol  scattering.  For  the  red  and  near-infrared 
region  of  the  spectrum,  variances  in  L(X)  are  proportional  to  changes  in  5a(^) 
and  therefore  aerosol  scattering.  Note  that  L(X)  in  the  channel  1  window  (short 
wavelengths,  small  particle  scattering)  is  greater  than  that  of  channel  2.  This  is 
primarily  due  to  the  slope  of  Lr(X)  and  leads  to  an  S 12  ratio  that  is  generally 
greater  than  1. 

B.  WATER  VAPOR  ABSORPTION  IN  THE  VISIBLE  /  NEAR 

INFRA-RED 

Although  aerosol  scattering  is  the  dominant  influence  on  satellite  retrieved 
radiances  at  visible  and  near-infrared  wavelengths,  the  role  of  atmospheric  water 
vapor  is  not  insignificant.  The  magnitudes  of  the  overall  L  in  these  wavelengths 
are  not  appreciably  effected,  but  often  the  goal  of  using  satellites  to  estimate 
aerosol  particle  characteristics  is  to  look  for  variations,  which  are  an  order  of 
magnitude  smaller  than  the  total  L. 

In  order  to  observe  sometimes  subtle,  yet  often  significant  trends  such  as 
global  climate  variations,  an  improvement  in  the  accuracy  of  aerosol  estimations 
is  necessary. 
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Fig.  2  NOAA-7  AVHRR  channel  1  and  2  response  functions  and  total 

atmospheric  radiance. 

Figure  3  shows  the  effect  of  water  vapor  alone  on  x  for  AVHRR  channels  1 
and  2.  Notice  that  for  channel  1,  x  is  greater  than  channel  2,  meaning  that  less 
absorption  is  taking  place.  This  would  tend  to  decrease  L2  more  than  Li  and 
thus  lead  to  an  over-estimate  of  the  S12.  Therefore  a  change  in  water  vapor  in  a 
region  may  be  misconstrued  as  a  change  in  aerosol  characteristics  or  disguise  a 
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significant  change.  This  implies  that  we  might  expect  to  see  a  minimum  Si2(w) 
at  w  =  zero,  and  an  increase  in  S12O)  with  w. 
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Fig.  3  NOAA-7  AVHRR  channel  1  and  2  response  functions  and 
atmospheric  transmittance  due  to  water  vapor. 

In  Figure  4,  the  influence  of  water  vapor  is  again  demonstrated  when  water 
vapor  is  added  to  LOWTRAN7’s  Tropical  atmosphere  model.  Notice  that 
while  radiances  decrease  in  both  channel  wavebands,  channel  2  is  more 
significantly  effected.  This  again  will  cause  S12  to  increase  with  w  even  though 
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Li  and  L2  decrease.  There  does  appear  to  be  some  increase  in  L  below  0.6  pm. 
This  is  most  likely  due  to  an  effective  shift  in  particle  size  due  to  coalescence  of 
water  vapor.  This  study  intends  to  look  at  the  influence  of  water  vapor  alone  on 


retrieved  radiances,  and  aerosol  growth  will  be  inhibited. 


Fig.  4  NOAA-7  AVHRR  channel  1  and  2  response  functions  and  the 
influence  of  water  vapor  on  atmospheric  radiance. 
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C.  BRIGHTNESS  TEMPERATURE  DIFFERENCE  (T4-T5) 

In  order  to  take  advantage  of  the  water  vapor  parameterization  of  the 
scattering  phase  function  proposed  in  this  study,  some  method  of  measuring  the 
water  vapor  present  would  be  required. 

A  significant  benefit  of  using  the  AVHRR  is  that  while  channels  1  and  2 
provide  information  on  optical  depth  variations,  channels  4  and  5  in  the  10  to  13 
|im  window  region  have  demonstrated  usefulness  in  atmospheric  water  vapor 
retrieval. 

Unlike  the  aerosol  backscattering  of  channels  1  and  2,  the  brightness 
temperature  is  the  result  of  single  direction  radiation  transfer,  from  the  sea 
surface  to  the  sensor.  Therefore,  x  is  directly  proportional  to  retrieved  T,  with 
absorption  reducing  the  radiance  reaching  the  sensor. 

From  Dalu,  et  al.  (1981): 

Ts  =  Ta  +  g(w)(Ta-Tp)  (2.8) 

Ts  =  Sea  Surface  Temperature  (SST) 

Ta  =  Brightness  Temp,  of  most  transparent  channel 
Tp  =  Brightness  Temp,  of  second  channel 
g(w)  =  function  of  water  vapor  content 
_ 1  -  Xq _ 

=  C  (1  -tp)  -(1  -  To) 

where  x«  and  xp  are  total  atmospheric  transmittances  for  channels  a  and  (5  and 
C=(TS- Tp)/(TS- Ta).  Solving  for  (Ta-T(j): 

(Ts*Ta) 

(Ta-Tp)  =  (2.9) 
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By  using  an  approximation  proposed  by  Prabhakara,  et  a/.(1974): 

(Ts'Ta)  =  (Ts-  T  a)Kw  sec0  (2.10) 

where  0  is  the  scanning  angle  and  K  is  an  absorption  coefficient.  Column  water 
vapor  (w  =  Jpwdz)  (g  m*2)  can  be  related  to  the  brightness  temperature 
difference: 


(TcrTp)  *  K  w  seed  (2.11) 

Although  both  g(w)  and  (Ts-  T  «)  increase  with  increased  water  vapor,  their 
ratio  remains  nearly  constant  (Dalu  1986).  This  leaves  us  with: 

w  =  A(Ta-Tp)cos0  (2.12) 

where  A  =  const.  =  (g(w)/(K(Ts-  T  a))- 

To  maximize  the  effectiveness  of  this  scheme.  Tot  should  be  minimally 

influenced  by  water  vapor  in  the  atmosphere  to  allow  the  retrieval  of  water 
vapor  temperatures  near  the  surface.  Tp,  on  the  other  hand,  should  be  more 

sensitive  to  the  atmospheric  water  vapor  distribution.  Under  such  condition:., 
(Tct'Tp)  would  indicate  the  influence  of  water  vapor,  which  would  be 

proportional  to  the  amount  (w). 

Channel  4  (~1  ljim)  and  channel  5  (~12|im)  windows  of  the  AVHRR  were 
designed  for  sensing  of  the  sea  surface  temperature  (SST)  corrected  for 
atmospheric  absorption  (Dalu  1986).  Despite  being  optimized  for  SST  retrieval, 
they  describe  a  linear  function  from  which  the  atmospheric  water  vapor 
correction  can  be  derived  (McMillin  and  Crosby  1984).  Figure  5  shows 
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NOAA-7  AVHRR  channels  4  and  5  response  function  curves  superimposed  on 
three  model  atmosphere  x  of  varying  water  vapor  amounts.  Notice  that  14  is 
generally  greater  than  T5,  leading  to  a  positive  T4-T5,  and  that  higher  T4-T5 


correlates  directly  with  higher  water  vapor  densities  (Tropical  >  US  Standard  > 
Midlatitude  Winter). 


Fig.  5  NOAA-7  AVHRR  channel  4  and  5  response  functions  and 
atmospheric  transmittance  due  to  water  vapor. 
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Although  not  the  ideal  sensor  windows  for  water  vapor  retrieval  (the 
channels  are  too  similar),  the  use  of  T4-T5  is  justified  to  within  the  accuracy  of 

the  approximations  made  to  achieve  Eqn.  2.12.  The  most  beneficial  aspect  of 
T4-T5  is  the  co-location  and  temporal  linkage  with  channels  1  and  2,  providing  a 

water  vapor  snapshot  simultaneously  and  using  the  same  path  as  the  aerosol 
optical  depth  (8)  measurements. 

Therefore,  an  estimate  of  vertically  integrated  water  vapor  density  (w)  is 
provide  by: 


w  =  A*(T4-T5)cos0  (2.13) 

Using  A*  =  19600  g  /(°K  m2)  and  a  scanning  angle  0  =  0,  Dalu  (1986)  applied 
this  relationship  to  a  radiative  transfer  model  using  a  variety  of  temperature  and 
humidity  profiles.  Figure  6  shows  the  correlation  coefficient  for  the  resulting 
line  is  R  =  0.99  with  the  error  given  as  ±  1.5  kg  nr2.  Comparing  results 
calculated  using  radiosonde  data  with  shipboard  measured  water  vapor  amounts. 
Figure  7  shows  the  correlation  coefficient  dropped  to  R  =  0.78  with  an  error  of 
±  4.0  kg  m*2.  The  decrease  in  accuracy  most  likely  resulting  from  calibration 
errors  as  well  as  inexact  temporal  and  spatial  correlation  of  ship  and  satellite 
data. 

At  channel  4  and  5  wavelengths,  scattering  is  a  relatively  insignificant 
influence  determining  brightness  temperature  (T).  The  dominant  factor  is  water 
vapor  absorption  of  surface  temperatures. 
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Fig.  6  Total  water  vapor  content  (kg/m2)  as  a  function  of  brightness 
temperature  difference  from  Dalu  (1986). 

Some  of  the  long  wave  radiation  emitted  by  the  sea  surface  is  absorbed  by 
water  vapor.  The  vapor  then  re-emits  radiation  based  on  temperature  in 
accordance  with  the  Planck  function.  If  the  atmosphere,  and  the  water  vapor 
included,  were  not  generally  cooler  than  the  sea  surface,  the  absorbed  radiation 
would  be  re-emitted  at  the  same  rate,  not  having  any  effect  on  radiance  or 
brightness  temperature. 
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Fig.  7  Comparison  of  water  vapor  content  from  ships  with  radiative 
transfer  model  runs  using  ship  radiosonde  data.  From  Dalu  (1986). 
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III.  PROCEDURE 


A.  INTRODUCTION 

This  study  uses  the  LOWTRAN7  model  to  improve  algorithms  used  in 
satellite  retrieved  optical  depth  measurements,  specifically  a  correction  due  to 
varying  amounts  of  water  vapor  present.  The  benefit  in  using  a  theoretical 
model  is  that  the  parameters  involved  can  be  specifically  quantified  and 
controlled.  Therefore,  a  precise  relationship  between  satellite  measured 
radiances  (L,  w/cm^-steradian-cnr1)  and  varied  parameters  such  as  optical  depth 
(8),  atmospheric  temperature  (T,  °K),  column  water  vapor  (w,  gm  nr 2)  can  be 
determined.  These  relationships  will  be  consistent  with  that  of  the  real 
atmosphere  to  the  degree  that  the  model  represents  the  real  atmosphere. 

B.  OVERVIEW 

Figure  8  is  a  flow  chart  of  the  procedure  beginning  with  a  user  defined 
atmosphere  of  known  quantities  of  aerosol  optical  depth  (8),  water  vapor 
densities  (pw),  and  surface  temperatures  (TSurf).  Three  programs  are  run  for 
each  data  set.  Both  LOWTRAN7  and  it's  filter  function  program,  LOWFIL7, 
make  use  of  formatted  files  (TAPEY.DAT)  for  input,  output,  and  to  pass  coded 
information  between  the  programs.  The  program  RETRIEVE  was  specifically 
designed  for  this  study,  to  compile  data  in  a  condensed  format. 

LOWTRAN7  is  supplied  a  modeled  atmosphere,  sun/satellite/earth 
geometry,  as  well  as  a  desired  wavenumber  (k)  band.  Wavenumber  dependent 
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radiances  are  calculated  considering  scattered,  absorbed,  and  emitted  radiation  as 
it  would  be  "seen"  by  a  perfect  sensor  at  the  satellite  location.  Figure  9  shows 
the  AVHRR  sensor  response  function  curves  (<D(k))  that  LOWFIL7  combines 
with  coded  information  derived  from  the  previous  LOWTRAN7  run  to 
calculate  the  effective  filtered  radiances  as  a  function  of  wavenumber  or 
wavelength.  The  sensor  weighted  total  integrated  radiance  is  then  determined. 

Ltot  =  |L(K)<I>(K)dK  (3.1) 

LOWFIL7  was  modified  for  use  in  this  study  to  calculate  the  integral  of  the 
sensor's  response  function  which  is  then  used  to  calculate  a  normalized  radiance 
(LOWFIL7A/B).2 


jLfltWtJdK 

L(chan)  =  f  (3-2) 

J<b(K)dK 

The  program  RETRIEVE  reads  the  output  files  from  LOWFIL7A/B  and 
writes  only  the  normalized  radiances  for  channels  1,  2,  4,  and  5  (Li,L2,L-4,L5). 
The  ratio  of  channel  1  to  channel  2  radiances  is  calculated  to  determine  the 
particle  size  index  parameter  (S12)  (Frost  1988).3 


o  _  (L)red  _  Li 
~  (L)nir  _  L2 


(3.3) 


2LOWFIL7  normalizes  by  ,  <t>max  instead  of  J<I>(K)dK. 

3Unlike  Frost,  this  study  does  not  first  remove  Lr. 
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Fig.  9  NOAA-7  AVHRR  niter  function  response  curves. 
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Channel  4  and  channel  5  radiances  are  converted  to  brightness  temperatures 

(Tb). 


Tb(4,5)  =  - CIk^ -  (3’4) 

+  L(chan  4,5)^ 

Where  Tb(4,5)  is  the  brightness  Temperature  (K)  for  the  radiance  L(4,5) 
(w/cm2-steradian-cm*l),  k  is  the  central  wavenumber  of  the  channel  filter 
(Kchan4  =  927.22  cm*1  and  Kchan5  =  840.872  cm*1),  and  Cl  and  C2  are  constants 
(Ci=  1.1910659xl0*12  w/cm2-steradian-cm*4  and  C2=  1.438833  cm*1  K) 

( NOAA  Polar  Orbiter  Data  1986).  Finally,  the  brightness  temperature 
difference  between  Tb4  and  TbS  is  calculated  (T4-T5). 

C.  INPUT  CHARACTERISTICS 

The  modeled  atmospheres  were  restricted  to  cloud-free  marine  environments 
where  backscattered  radiance  is  dominated  by  aerosol  particles  (Durkee  1984). 
Since  the  satellite  senses  the  effect  of  vertically  integrated  parameters,  the  models 
used  were  homogeneous  from  the  surface  to  the  tropopause.  For  ease  of 
integration,  the  model  tropopause  was  set  to  10  kilometers  (the  average  height  of 
the  mid-latitude  tropopause  is  about  11  kilometers  (Stull  1988)).  From  the 
tropopause  to  the  satellite  height,  only  background  extinction  due  to  meteoric 
dust  is  represented  from  30  to  100  kilometers.  The  following  are  the  important 


user-defined  parameters: 

aerosol  optical  depth  (8a(X))  . varied 

column  water  vapor  (w)  . varied 

tropospheric  air  temperature  (Tair)  . 260  K 

surface  temperature  (Tsurf)  . varied 


24 


surface  albedo  (Salb)  . 0.0  (blackbody) 

satellite  height  (H)  . 850  km 

season  . spring-summer 

marine  phase  function  (P(0)) . from  Shettle  and  Fenn  (1979) 

wavenumber  windows,  resolution 

chan  1,2  . 8400-20050  cm-1,  50  cm*1 

chan  4,5  . 770-1005  cm*1,  5  cm*1 

Molecular  constituents  such  as  O3,  CH4,  N2O,  CO,  CO2,  O2,  NO,  SO2,  NO2, 
NH3,  and  HNO3,  used  LOWTRAN7  model  default  profiles  for  midlatitude 


summer  and  were  kept  constant  for  all  model  runs.  The  sun/satellite  geometry 
remained  the  same  throughout  all  the  model  runs  with  the  satellite  and  sun  in  a 
straight  line  perpendicular  to  the  earth's  surface. 

1.  Optical  Depth  Profiles 

Satellite  windows  in  the  red-visible  and  near-infrared  sense  the  result  of 
vertically  integrated  parameters.  The  most  significant  influence  on  retrieved 
radiances  is  the  vertical  integral  of  aext  or  8.  To  study  this  effect,  ten  profiles 
of  8  with  constant  oext  to  10  km  were  designed. 

Using  the  relationship 


(Jscat06^'^  (3.5) 

and  an  oScat  of  .05  km*1  at  X  -  .55  pm,  a  profile  of  Gscat  as  a  function  of  X  was 
derived  (Oscat(^))*  After  integrating  over  the  10  km  height  to  determine  8, 
Figure  10  shows  this  curve  (M2)  and  four  multiples  of  this  profile  that  complete 
the  first  set  of  optical  depth  profiles. 
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Fig.  10  Multiple  optical  depth  profiles  (M). 

Figure  1 1  shows  a  second  set  of  profiles  created  using  the  same  relationship  and 
adding  equal  quantities  to  create  parallel  profiles. 

2 .  Water  Vapor 

All  the  water  vapor  was  contained  in  the  homogeneous  troposphere  (0 
to  10  km).  The  maximum  value  was  selected  from  the  tropical  atmosphere 
model  included  in  LOWTRAN7  (3.322  gm  cm*2).  User-defined  profiles  are 
input  as  water  vapor  mass  density  (pw)  (gm  m'3),  which  after  integrating  in  the 
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Fig.  11  Parallel  optical  depth  profiles  (P). 
vertical  gives  column  water  vapor  (w)  of  equivalent  magnitudes  in  units  of  gm 

cm*2. 

H 

w  =  Jpw  dz  (3.6) 

Ten  intermediate  quantities  between  0  and  3.322  gm  cm'2  were  selected  with  0 
included  as  the  eleventh. 
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3.  Surface  Temperature 

T4-T5  is  a  measure  of  water  vapor  (Jacobowitz  and  Gruber  1990)  and 
will  be  compared  with  the  user  input  water  vapor  concentrations,  as  an  empirical 
method  of  estimating  w.  T5  is  more  sensitive  to  water  vapor  effects  on 
atmospheric  transmittance  while  T4  senses  more  of  the  near-blackbody,  long* 
wave  radiation  from  the  surface.  In  order  to  demonstrate  relationships,  the 
model  required  a  warmer  surface  temperature  than  the  atmosphere  when  the 
troposphere  is  isothermal.  With  equal  atmospheric  temperature  (Tatm)  and 
surface  temperature  (Tsurf),  the  measured  brightness  temperatures  would  be 
identical.  Also,  surface  temperatures  colder  than  the  temperature  at  10  km  are 
unrealistic.  Therefore,  while  the  tropospheric  temperature  was  held  constant, 
the  surface  temperature  was  varied  to  achieve  four  representative  modeled 
Tatm/Tsurf  relationships: 

Tatm . 260°K 

Tsurf(w) . 270°K 

TsurfOO . 280°K 

TSUrf<y) . 290°K 

TSUrf(z) . 300°K 

D.  PROCEDURAL  SUMMARY 

The  model  was  run  at  the  Naval  Postgraduate  School’s  (NPS)  Interactive 
Environmental  Digital  Analysis  Laboratory  (IDEA  LAB).  A  complete  set  of 
input  files  consisted  of  two  sets  of  five  8  profiles,  with  eleven  w  variations,  four 
surface  temperature  conditions,  and  two  wavenumber  windows.  These  880  files 
(2X5X  11  X4X2)  generated  1760  intermediate  and  final  output  files,  with 
the  significant  data  ’retrieved’  and  collated  in  one  resultant  array.  All  the 
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calculations  were  done  as  a  large  batch  job  requiring  about  4.5  hours  to  run. 
The  job  was  completed  in  both  the  single  and  multiple  aerosol  scattering  modes 
of  LOWTRAN7. 
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IV.  RESULTS 


A.  THE  EFFECT  OF  AEROSOL  OPTICAL  DEPTH  VARIATIONS 
ON  S12 

In  general,  an  increase  in  wavelength  dependant  aerosol  optical  depth 
(Sa(^.))  resulted  in  an  increase  in  Li  and  L2  while  S12  values  decreased.  Since 
S12  is  generally  greater  than  one  (Li>L2).  Increasing  the  radiance  in  both 
sensor  windows  by  similar  magnitudes  has  a  greater  influence  on  the  smaller  L, 
driving  S12  toward  one.  This  relationship  will  remain  true  because  of  the 
Rayleigh  contribution,  which  is  greater  for  Lj. 

1.  Multiple  Sa(A.)  Scheme 

Increasing  5a  profiles  result  in  decreasing  S12  ratios,  with  the  rate  of 
decrease  diminishing  due  to  the  effect  of  increasing  8a(X)  slopes.  Figure  12 
shows  the  S12  variations  resulting  from  the  multiple  5a  profiles  of  Fig.  10. 

2.  Parallel  5a(X)  Scheme 

By  constraining  the  curves  to  be  parallel,  an  increase  in  5a  is  an  equal 
increase  at  all  wavelengths.  The  S12  ratios  decrease  at :  ore  rapid  rate  than  in 
the  previous  scheme  since  a  change  in  slope  is  no  longer  limiting  the  decrease  in 
Sl2.  Figure  13  shows  the  S12  variations  resulting  from  the  parallel  5a  profiles 
of  Fig.  11. 

B.  THE  EFFECT  OF  SURFACE  TEMPERATURE  VARIATIONS  ON 
THE  T4-T5  REPRESENTATION  OF  W 

Although  the  modeling  accomplished  in  this  study  was  not  designed  to  prove 
Dalu's  technique  for  determining  vertically  integrated  water  vapor  density  (w) 
from  T4-T5,  some  significant  aspects  of  his  theory  are  supported  by  the  results. 
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COLUMN  WATER  vapor  ((cm -2) 
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Fig.  14  T4-T5  versus  column  water  vapor  (w)  for  a  variety  of 

surface  temperatures. 
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TABLE  2 

MODEL  OUTPUT  COEFFICIENTS  APPLIED  TO  DALU'S 
THEORY  (w  =  A(T4-T5)) 


Tair/Tsurf  (K/K) 

A  (g  K*1  nr1) 

260/300 

.7615 

260/290 

.9944 

260/280 

1.4995 

260/270 

3.6653 

Dalu 

1.9600 

C.  THE  W  CORRECTION  TO  S n 

In  order  to  present  the  correction  to  satellite  retrieved  S12  values  required 
due  to  the  presence  of  water  vapor,  the  results  will  be  displayed  first  graphically, 
then  in  a  new  parameterization  that  can  be  applied  to  the  aerosol  scattering  phase 
functions  used  in  retrieval  of  aerosol  characteristics  from  satellites. 

1 .  Graphical  Application 

RETRIEVE  lists  the  following  results  for  each  modeled  atmosphere: 

Input  Parameters:  Output  Variabl££l 

8a(X)  L2 }  =>Sl2 

Tsurf 

W  Tbsl  =>T4'TS 

A  set  of  output  variables  was  calculated  for  each  variation  of  input  parameters. 
The  combined  results  are  displayed  in  Fig.  15. 
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To  use  this  figure  to  determine  S12  corrected  for  w:  Begin  with  NOAA-7  AVHRR  retrieved  data  (S12, 
74*75).  Read  74-T5  value  on  lower  half  of  plot  and  look  for  the  intersection  with  Dalu  curve.  This 
will  determine  column  water  vapor.  Move  up  (parallel  to  the  grid  lines),  to  the  satellite  measured  S12 
value.  Follow  the  contours  of  the  nearest  curve  to  the  Sj2  axis  and  read  Si2(dry).  This  will  be  the  S 12 
with  the  effects  of  water  vapor  in  the  column  removed. 


0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5 

Column  Water  Vapor  (g  cm"2) 


Fig.  15  Si2  correction  for  column  water  vapor. 
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Both  the  multiple  and  parallel  S12  results  are  displayed  and  seem  to  follow  a 
similar  relationship  to  w.  Note  that  M2  and  PI  curves  are  coincident.  They  use 
the  same  input  5a(X)  values  that  formed  the  basis  for  each  set  of  curves. 

2.  The  Reformulated  Parameterization 

Looking  at  the  S12  values  plotted  for  varying  w  and  5a(A.)  curves, 
there  appears  to  be  a  functional  relationship.  It  is  apparent  that  the  water  vapor 
influence  of  S12  is  not  dependent  on  the  shape  of  the  optical  depth  profile,  but 
was  a  function  of  w  and  the  magnitude  of  S12  when  w  =  0  (Si2(dry)).  Figure 
16  shows  that  when  the  ratio  of  Si2:Si2(dry)  was  plotted  as  a  function  of  w,  the 
shape  of  the  curve  was  conserved  for  all  model  outputs  (multiple  and  parallel). 
The  slight  variation  at  higher  w  is  partly  due  to  computational  round  off  errors. 
The  curve  that  best  fits  the  data  is: 

Si2(dry)  =  (1  +  .0332  Vw)  S12  (4.1) 

with  r^  =  0.995. 

From  Eqn.  4.1,  measured  S12  and  estimated  w  from  T4-T5  determines 
Si2(dry). 

D.  SINGLE  VERSUS  MULTIPLE  SCATTERING 

All  data  sets  were  run  with  both  single  scatter  and  multiple  scatter 
approximations.  The  results  were  calculated  to  the  seventh  decimal  place  and 
were  precisely  the  same  for  either  case.  The  model,  therefore,  supports  Frost’s 
use  of  a  simplified  form  of  the  Radiative  Transfer  Equation  (RTE)  that  assumes 
single  scatter  (Eqn.  2.1). 
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V.  CONCLUSIONS 


Present  techniques  for  satellite  retrieval  of  aerosol  characteristics  using  the 
AVHRR  carried  on  the  NOAA  series  platforms  require  that  analyzed  pixels  be 
free  of  clouds.  This  is  reasonable  since  liquid  water  effects  on  radiative 
processes  are  quite  variable  and  difficult  to  measure.  Clouds  at  any  level  in  a 
satellite  scene  can  obscure  the  sources  it  was  designed  to  sense,  or  saturate  the 
sensor  with  reflected  radiance.  It  is  understood  that  water  in  a  gaseous  form 
would  still  be  present,  but  its  influence  on  L  in  the  wavebands  most  effectively 
used  to  estimate  aerosol  parameters  (Chan.  1  and  Chan.  2)  has  been  assumed  to 
be  insignificant  in  previous  studies. 

Although  aerosol  scattering  is  still  the  dominant  influence  on  L  in  these 
wavebands,  this  study  has  shown  that  water  vapor  does  introduce  an  error,  and 
that  the  error  can  be  removed.  This  does  not  invalidate  the  present  techniques, 
but  provides  a  significant  improvement  that  can  be  applied  through  the  use  of  the 
variable  scattering  phase  function  introduced  by  Frost  (1988).  Since  aerosol 
variations  can  be  of  global  scale,  even  very  small  changes  may  have  global 
ramifications. 

Through  the  use  of  the  LOWTRAN7  atmospheric  propagation  model,  this 
study  has  shown  support  of  the  following  techniques  and  theories  presently 
accepted: 

•  The  concepts  of  Frost's  particle  size  index  parameter  (S12)  using 
AVHRR  channels  1  and  2. 

•  Frost's  use  of  a  simplified  radiative  transfer  equation  based  upon  an 
assumption  of  single  scatter. 

•  Dalu's  water  vapor  retrieval  technique  using  AVHRR  channels  4  and  5 
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By  introducing  varying  quantities  of  water  vapor  to  the  modeled 
atmospheres  of  differing  optical  depth  (S)  and  Tsurf/Tajr  profiles,  additional 
information  has  been  generated: 

•  The  error  introduced  to  S12  due  to  water  vapor  is  demonstrated  and 
quantified. 

•  The  error  has  been  shown  to  be  a  function  of  Si2(dry)  and  column 
water  vapor  (w). 

The  influence  of  optical  depth  variations  can  be  isolated  from  the  variation 
due  to  w  through  the  use  of  the  variable  scattering  phase  function  (P(©)) 
proposed  by  Frost,  (1988).  Using  Eqn.  2.1,  S12  with  2.0gm  cnr^  of  water 
vapor  produces  a  15%  error  in  P(0)  and  therefore  retrieved  aerosol  optical 
depth  (8a). 

By  using  a  model  to  determine  water  vapor’s  influence  on  satellite-measured 
aerosol  characteristics,  the  parameters  normally  estimated  can  be  strictly 
controlled.  As  known  quantities,  their  theoretical  relationships  to  satellite- 
measured  quantities  can  be  observed.  The  next  step  is  to  apply  the  correction 
discovered  here  to  empirical  satellite  data  with  verifying  in  situ  data  sources  such 
as  aircraft  or  ships. 
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